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Abstract 

In this paper we employ the continuum approximation of Dyson to determine the 
asymptotic gap formation probability in the spectrum of N x N Hermitean matrices 
associated with orthogonal polynomials that are solutions of indeterminate moment 
problems. 
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I Introduction 

Random matrix ensembles originally proposed by Wigner as a phenomenological descrip- 
tion of the statistical properties of the energy spectra fl]],@ has recently seen applications 
in other area of physics such as Quantum Chaos 0, transport in disordered systems Q and 
2-d quantum gravity [f|. From the random matrix point of view|J,@, the quantities of 
interest (in increasing order of refinement) are the eigenvalue density, the density-density 
correlation function and Ep[n, J] — the probability that an interval, J, contains exactly n 
eigenvalues. [This interval can be generalized to the case where it is the union of several 
mutually disjoint intervals f§.||]. 

This paper shall focus on Ep[Q, J], henceforth denoted as Ep[J\. The probability that there 
is no eigenvalue in a scaled interval (—t,t), is kiEp(— t, t) ~ —^fit 2 — (1 — (3/2)t; first 
determined by Dyson|[U|, using the methods of classical electrostatics, potential theory 
and thermodynamics, for random ensembles with unitary (j3 = 2), orthogonal {(3 = 1) 
and symplectic (/3 = 4) symmetries. The missing constant and O(lnt) terms were later 



found in 0,0 and |3| 



For Hermitean matrices (j3 = 2) a fundamental result of Gaudin and Mehta shows that 

E 2 [J]=det(l-Kj), (1) 
where Kj is the kernel of an integral operator K(x, y) over J; 

K(x, y) = e-^VV^/ 2 Y, Pn (x)Pn(y), (2) 

n=0 

dx w{x)p n (x)p m {x) = 5 mn , w(x) := e~ u{x) , (3) 

where P is the interval available to the eigenvalues. Here u(x) in the Coulomb gas 
interpretation is an external (confining) potential that holds together a gas of repelling 



particles. An alternative representation for the Fredholm determinant is 



recognized as the ratio of two partition functions, where 

(5) 



and 

N 

W = -P ln|x a -x 6 | + ^w(x a ). (6) 

0<a<b<N a=0 

A Fredholm determinant of this kind also appears in context of the physics of interacting 
bosons in one dimension ]14] where it gives the probability of emptyness formation of a 



segment of the real line. 

For sufficiently large N it is expected that the discrete Coulomb gas can be well approx- 
imated by a continuum Coulomb Fluid[|]l| with an equilibrium density cr(x) of charged 
particles that satisfies a Hiickel-like self consistent equation; 



u(x) — (3 J dya(y) In \x — y\ + ^1 — In [cr(x)] = A = constant, 



(7) 



and gives the free energy, at equilibrium, with precisely N charges (eigenvalues) contained 
in an interval J; 



F[J] = ]-AN + i f dxu(x)a(x) - \ ( 1 - ^ ] f dxa(x) In 



a(x)} 



subject to Jjdxa(x) = N. This approach has the advantage of being more general and 
has been shown to supply accurate asymptotics of the level-spacing distribution for the 



Laguerre ensemble at the hard edge[16 for all (3 > 0. We expect that the continuum 
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Coulomb Fluid description will also be useful in determining Ep[J] of orthogonal polyno- 
mial random matrix ensembles associated with non-classical weight, w(x) = e~ u ^ x \ where 
a explicit form of the polynomials are not known. 

This paper is organized as follows: In section II, we shall introduce and discuss the 
properties of Hermitian random ensembles (hence forth called the a— ensembles) with 
u(x) = x a /a,( a > 0, x G (0, oo)). It will be shown in the continuum approach the 
parameter a plays the role of a coupling constant with a critical value (oe c ), above which 



the density is universal (i.e. a independent) at both the hard and soft edges [[L8[]. For 
a > a c , the density at the soft edge remains universal, however, an explicit a dependence 
is found at the hard edge. 



By using a heuristic argument based on the leading asymptotics found by Dyson (JTUJ fol- 
low by an un-folding transformation of the density, we obtained an asymptotic formula 
for, E 2 [J], the probability that the interval J contains no eigenvalue. This is first used 



to test against known results |[18fl , ||16[ and later make predictions on the gap formation 
probability of the a ensembles. 

In section III a screening theory of the Coulomb Fluid is presented to justify the results 
obtained in section II. Level Spacing Distribution for a = 1/2 is case presented in detail 
in section IV. This paper concludes with section V, in which we give E 2 [J] for certain 
random matrix ensemble arising from the double scaling limit of two-dimensional quan- 
tum gravity. 

II Behaviour of density at edges 
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As we shall be interested exclusively in ft = 2 case, Eq.(7) becomes linear and reads; 



u(x) —ft dya(y) In \x — y\ = A, 
Jo 



(9) 



where b, the upper band edge, is related to the total number of particle (i.e. eigenvalues) 
through Jq dya(y) = N. By taking a derivative with respect to x, Eq.(9) is converted into 
an integral equation with Cauchy kernel; 



u'(x) = ftP t dy 
Jo 



x-y'' 



(10) 



whose solution can be obtained by a standard technique p9| ; 

, , 1 



a(x 



ft^ 





f b dy 


/ y 




o y — x y 


b-y 



u'{y). 



'IV 



Since 



dy 



(12) 



lo X- V> Jy(b-y) 

the general solution of Eq.(10) is found by adding to Eq.(ll) the solution of the homoge- 
neous equation, c ° nstant . However, it is clear that by including the homogeneous solution, 

y/x{b-x) 

the equilibrium free energy will be larger. The homogeneous solution will not be included. 
It can also be verified by direct computation that the density obtained without the homo- 
geneous solution compares well with the density obtained from orthogonal polynomials. 
A generalization to the case of finitely many non-overlapping has been studied 



by Akhiezer [19 



u'(x) = ftP I dy^L, xeJ 



j x — y 



(13) 



where 



J — Uj! =1 Jfc, Jk — (dk,bk), ak < bk- 



(14) 



The solution to Eq.(13) is 

a(x) 



7T 2 ' 



fl( bj —^]u'(y), xeJ. (15) 



2 13 ^ ^ \x- a k ) Jj y - x\^fJi\y - a t 
This integral representation can be useful in determining the level spacing distribution in 
the multiple interval case. 
If we now take u(x) = - x a , a > then 



b — x 


f b dy 


/ y 


/ x J 


o y — x y 


b-y 



*w = -tz\i — i y a 1 



b 



a-l 



2r(l/2)r(a + 1/2) b-x „, 
^ V [a) i— m- a ,l,3/2,l-x/6), (16) 

where F(a, 6, c, x) = 2-^1 (°> c i s ) is the Hypergeometric function. Note that for a — 1, 



the result for the Laguerre Ensemble is recovered; a[x) = ^y^jr- We have in the Hy- 
pergeometric function the following parameters, a = 1 — a, b = 1, c = 3/2, and a + 
b — c = ~ — a. From this the behaviour of the density at the the hard (x ~ 0) and soft 
(x ~ 6) edges fl8| can be determined. Consider a > 1/2, and 1 C &, we find from the 



Gauss summation formula, 

r(3/2)r(a - 1/2) 



F(l- a, 1,3/2,1) 



r(l/2)I> + l/2)' 



^-iW^P^f?'** (17) 

The density has a square root singularity at x = 0, independent of a and is therefore 
universal. Near the soft edge, x ~ b, 



, x 1 r(a+l/2) in , , b-x 1 

which shows that the density is universal near b. However, for \a\ < 1/2, a completely 
different behaviour at the hard edge occurs, this can be seen as follows: A translation,^ — > 
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1 — x gives 



F(i - «, i, 3/2, i - X ) = rja + i/^rci/^ ^ " a ' 11 3/2 " a ' x) 
+ ^-i/2 r (3/2)r(i/2-c) F + q + 

1 (1 — a) 

In view of the combination of the parameters (a, 6, c), the second term of the right hand 
side of Eq.(19) dominates as x — > 0. Putting x = in the second Hypergeometric function; 
we find F(a, b, c, 0) = 1, < a + b — c < 1. Therefore, 

. , tan7ra 1 . , 

~ ^> X « & - (20) 

Eq.(20) has two features that are worth noting; 1. a(x) near the hard edge is independent 
of the macroscopic parameter b [which increases with N]. This type of behaviour at 
the hard edge is also be observed in the unitary q— Laguerre Ensemble with u(x; q) = 
En=oW + (1 - q)xq n ], < q < 1 x > 0; 

l-q N 



a N (Q) 



which rapidly becomes N independent as increases ||20|| , 2. The hard edge density has 
a explicit a dependence. 

We remark here that the kernel in the unitary q— Laguerre Ensemble reads, 

K(x,y) 



constant 
x-y 



' \ V4 / \ / \ / x 1/4 , 



sin — In a; cos — In y — — sin — In y } cos — In x 
V V 2 7 J Vl J l2 7 y ) 12/3 , 



(21) 



in which 7 := ln(l/g) ^> 1 and — > 00 [21]. Therefore the so-called universal density- 



density correlation function at the hard edge < p(x)p(y) > — < p(x) >< p(y) >= S(x — 



y) < p{x) > —[K(x, y)} 2 = (xy)~ 1/2 (x 1/2 + y 1/2 )/(x - y) 2 @ is seen to be violated. This 
is due to the weakly confining nature of the q— Laguerre potential; u(x; q) ~ [In a;] 2 , x —>■ 
oo, producing a heavily depleted density at the origin. We expect in the a— ensemble a 
similar phenomena to occur for a < 1/2. In principle, the gap formation probability of 
the unitary a— ensemble at the hard edge can be computed directly from the Fredholm 
determinant in Eq.(l), however, the associated orthogonal polynomials being related to 



an indeterminate moment problem[25. 23] is not known, unlike the unitary q— Laguerre 



Ensemble [24]] . A direct evaluation of the Fredholm determinant would be difficult. To 
extract the leading term in — In E 2 (0, s) we shall adopt the Coulomb Fluid approximation 
mentioned in section I. 

Before we give results on the a— ensemble, it is interesting to note that for a = 1/2, 



<T(X) 



7T 2 /3 y/x 



In 



1 + yJl-x/b 

1 - Jl - x/b 



b = 7i 2 N 2 



(22) 



where we have made use of the identity 

F(l/2,l,3/2,x) 

to arrive at Eq.(22). 
Thus for x <C b 



2^ 



In 



l + y/x 
1- 



a(x) 



1 ln(4b/x) 



7T 2 f3 \fx 



(23) 



while for x ~ b 



a(x) ~ 



n 2 (3 



(l_x/6) 1 /2 + £(l-x/6)S/a + - 



remains to be universal. We may therefore interpret a = 1/2 to be the critical point at 
which there is a logarithmic correction to the density. 
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In order to gain insight into the possible leading terms in the gap formation proba- 



bility, we examine the results first obtained by Dyson for the Circular ensemble JTO 
where by construction the density is uniform over the circle, N/2n. It was shown that 
— lnE 2 (0, s) ~ constant s 2 , s > 1. Here s being expressed in terms of the average 
distance between the charges is dimensionless. The leading term is simply proportional 
to the square of the number of particles excluded in an interval of length s. Using the 
unfolding transformation in the density, we can always define a new variable y, such that 
a(x)dx = dy, or equivalently ^ = a(x). The density in y is unity. At this point, we may 
make use of Dyson's asymptotic and conclude that the leading term is proportional to 
the square of the number of particles excluded in that interval. The probability that the 
interval J contains no eigenvalues is 



j 



dxa(x) 



2 



(24) 



-\nE p [J}~ 

Based on Eq.(24), the probability that an interval of length a (measured from the origin 
contains no eigenvalue is, 

^H-e-^ 2 , 

where 

J\f(a) : = / dxa(x), 
Jo 



is the number of charges excluded in (0, a). Taking a(x) = A y ^r, we find 



b f a / b 1 1 — x 

dxa(x) 



, 1 (abfl 2 



(25) 



lift JO V X 7T/3 

In the limit b — ► oo (or iV — > oo, since b oc N) and a — > such that the combination 
ab = s is finite, gives 

-In £7/9(0,5) ~ s, (26) 
9 



agrees with the know leading behaviour [|I8|1 , |26| , fl6|| . By taking J = (a, 6), a «b and 
scaling b — a with respect to the true soft edge density ctn(N) ~ iV -1 / 3 ~ b^ 1 ^ 3 , we find, 

-ln^(0, S )~ S 3 , where s := , (27) 

which agrees with that found in[p^|,||27|. Observe that although the density in the con- 
tinuum approach vanishes at b, the correct dependence at the soft edge, i.e., a(b) ~ 



is found [[26 



For the a— ensembles, 

-ln£^ a) (0,a) ~ a 2a , < a < 1/2, (28) 

while a = 1/2, 

- In 4" 2 » (0, a) ~ f" dJ^H ~ h, f ^ V (29) 
p Jo \/x \ a 



II Screening theory 

In this section we shall consider the unitary case {(3 = 2). With a slight change of notation, 
the free energy of N particles restricted in an interval / = (0, b) under the influence of an 
external potential u{x)/ (3 is 

F[a,I,N] = J dxu(x)a(x,I,N) - J^dx J dya{x,I,N)\n\x -y\a(y,I,N), (30) 

where a(x, I, N) is the equilibrium density that minimizes F[a, I, N], subject to 

a(x)dx = N. 



i 

According to Eq.(4) the quantity of interest is the change in free energy, i.e. the free 
energy when the charges are excluded from the sub- interval J = (0, a), a < b of I minus 
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the free energy when the charges are in /: 



AF[J, N) = F[J C , N] - F[I, N], 



(31) 



in the limit N — > oo. Here the interval J c is the the interval where the charges are 
distributed, when there are no charges in J. For later use we introduce, 



Aa{x, J, N) := a(x, J c , N) - a(x, I, N), 



(32) 



and require that all the cr's are zero outside their respective supports. The support of 
Aa(x, J, N) is L = I U J c . 

Using the definitions of AF and Act in Eqs.(30) and (31), a simple calculation gives, 



AF[J,N] = J L dx u(x) - 2 £dy\n\x - y\a(y, I , N) 



Aa(x, J, N) 



- f dx f dyAa(x, J, N) In \x - y\Aa(y, J, N). 

2 J L J L 



(33) 



Since a(x, I, N) is the density that minimizes the free energy; [•••], in Eq.(33) is equal to 
the chemical potential, /i(J, N), for x E I and is a function of of x, for x E L — I = J c — /. 
Charge neutrality states that J L dx Aa(x, J, N) = 0. Therefore the first term in Eq.(33) 
can be rewritten as 



J dx u(x) — 2 J dy In \x — y\cr(y : I, N) 



Aa(x, J, N) 



dx 



u 



(x)-2j i dy In \x - y\a(y, /, N) - /i(J, iV)j Aa(x, J, N) 

= f dx u(x)-2 f dy\n\x - y\a(y,I,N) - /j,(I,N) Aa(x, J, N). 

The term [• ■ •] under J L _j dx is of the same order of magnitude as the change in the 
chemical potential when the particles are excluded from J: A/x( J, N) = fi(J c , N) — 
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p(I, N). A/i( J, iV) decreases as N grows and in the limit, iV — > oo, it is proportional to the 
number of particles excluded from J. Therefore the first term in Eq.(33) is approximately 
Ap( J, N) times the number of particles in J c — I, which is the number of particles that 
spilled over when all particles are excluded from J. In the bulk and at the hard edge 
this number goes to zero when N — > oo and therefore the first term in Eq.(33) can be 
neglected. At the soft edge the first term contributes a term quadratic in the number of 
particles excluded from J. This is of the expected form from the unfolding transformation 
and will not affect our final result. We therefore have, 



here Aa(x, J, N) is simply —cr(x, I, N) for x G J. For x £ J, the expression is more 
complicated. From now on we shall work exclusively with rescaled quantities such that 
lengths are measured in the units of \J\(— a). We also translate the origin of the real axis 
to the left end of J. With the scaling, x — > ax, and y — > ay, and the introduction of p(x) 
and X(x), 




dx dyAa(x, J, N) In \x — y\ Aa(y, J, N), 



(34) 



p(x) :=p(x,J,N) : 



aa(ax), x G (0, 1) 



(35) 



\{x) := X(x, J, N) = aA(ax, J,N), x G [1, b/a). 



(36) 



The change in free energy when expressed in terms of scaled quantities reads, 





dx i dyAa(ax, a) In \ax — ay\ Aa {ay, a) 



— / dxAa(x, J) — — / dx dyAa(ax, a) In \x — y\ A(ay, a) 

2 \_J L J 2 J LI a JL/a 




\ f dx f dy\(x)\n\x - y\X{y) - J / dx f dyp(x) In \x - y\p(y). (37) 

2 Jjcja JJ c /a 2 JO JO 
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Note that the first integral in the second equality of Eq. (37) vanishes due to charge 
neutrality,/^ dxAa(x, J) = 0. The problem of finding Aa(x, J, N) in the limit iV — > oo 
can be reformulated in physical terms. The interval I is force (or field) free since J is a 
conducting region. Outside / we have a potential that grows when / tends to infinity. We 
introduce negative charges in J with density — a(x, I, N) and precisely the same numbers 
of positive charges outside J. [The is the charge neutrality condition stated above.] The 
charges redistribute themselves in such a way that the interval L — J is force free again. 
The charges in J are screened by those outside J. This screening process is not as efficient 
as in three dimension since the charges are only allowed to be built up along a line while 
the force extends out in a plane. X(x, J) is then determined by the condition that the 
total force in J c /a is zero. 
The force produced by the charges in J is 

f(x)=[ 1 dy-^. (38) 
Jo x — y 

The force produced by the charges in J c is 

Mx) = f dy^. (39) 

J J c /a x — y 

Therefore the total force is 

fjc(x) + f(x) + f 1 (x) = 0, xeJ c /a, (40) 

where // is the combined force in J, produced by the external potential u{x)/2 and the 
charge distribution Aa(x, I, N). 

Two constraints are required to solve Eq.(40): 1. The number of charges in J c /a equals 
that in J/ a, 2. Aa(x, J c ) is positive semi-definite for x G J c , and minimizes the free energy. 
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We shall also require that at an end point of J c / a the density vanishes. The reason is 
simply that the charges would otherwise distribute in a larger interval. This argument 
can be generalized to the case where J c = lf k=1 Ik, with the additional requirement that 
the chemical potentials to be the same in each intervals l k . The above formulation can 
be applied to X(x). In the above scheme we make the following approximation: If J is 
not at the end of / where new charges gather, the total potential can be approximated 
by a constant all the way up to infinity and not only inside /. Observe that in this 
approximation A is supported everywhere on the real axis outside J; except in the hard 
edge case where it is only a half line. Call this interval P. In the rescaled units we 
have as an example in the hard edge problem, P = (l,oo), while in the bulk problem, 
P = (—00, 0) U (1, 00). In terms of the interval P, the change in free energy reads, 



To proceed further, we can express A in terms of p by the use of the condition that the 
total force is in P. As above we find that the total force is 



The solution of Eq.(42) depends only on P and \\ x edP- For, P = (1, 00), the solution is 



The most important feature about this solution is that A is a linear functional of / and 
therefore also of p, (see Eq.(40)), which does not depend on a and N. To proceed further, 




^ J dx J dy\(x) hi \x-y\\(x) - dx dyp(x) In \x - y\p(y). (41) 




(42) 




(43) 
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we now develop / in a "multipole" (or moment) expansion; 

/(*)= f'dy^-, xe(l,oo) 
Jo x — y 

oo i ! 

= E I dvv m p(y)- m 

m=0 X JU 

The moments 

A rn := [ dxx m p(x), (45) 
Jo 

will play an important role in a later development. Since A is a linear functional of / and 
/ is a linear functional of A' m s, we conclude that A is a linear functional of the A' m s. By 
expanding y/xp(x), x G [0, 1], in terms of Legendre polynomials P n (2x — 1), 

oo 

Vx~p(x) = ]T r m P m (2x - 1), x £ [0, 1] (46) 

m=0 

we find also that ^p(x) is a linear functional of the moments A m+ i. In summary, the 
densities can be expressed as linear functional of A m and A m+ i. Since the change in the 
free energy is quadratic in the densities, it can be expressed as a sum quadratic in the 
A' m s. A simple estimate gives, 

|An| < IM (47) 

from which we deduce that the change in free energy is bounded from above by A%. This 
we recognized as the square of the number of charges excluded in (0, a). Precisely the 
same analysis can be carried through in the bulk scaling case. 

The upshot of the above analysis is that AF is quadratic in the moments of the local 
densities, when we adopt the Coulomb fluid approximation. From this we may infer two 
consequences: 1. The probability of gap formation depends only on the local densities. In 
the bulk scaling case the density is always a constant and the asymptotics level spacing 
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distribution found by Dyson is reproduced, 2. Higher moments are nominally linear 
functions of the number of excluded charges [a logarithmic dependence is in general missed 
in the continuum approach, because AF(J) disappears when the length of J tends to 
zero.]. Therefore, AF for large spacing is proportional to the square of the number of 
excluded charges in the appropriate interval. On the other hand it is shown thatjl6| the 



gap formation distribution of the Laguerre Ensemble at the hard edge, the large interval 
formula involves quadratic, linear and a logarithm of the number of excluded charges 
(oc s/s). The breakdown of the present theory can be traced back to the fact that the 
density of the Laguerre Ensemble contains a point charge at the origin and therefore can 
not be adequately described by its moments. We may be conclude that leading term of 
— hi Ep[ J], is proportional to Af 2 [J], where M[J] is the number of charges excluded in J. 
This justifies the un- folding transformation discussed in section II. 

We now proceed to express the change in free energy in terms of A' m s; substitute Eq.(43) 
into Eq.(42), we have 

u \ 1 1 r d v i — 7^ 2A ™ 



7T 2 y/x-1 Ji x-y v ^ Q (2y-1 



,m+l 



g 2A m 5(1/2, m + 1/2) F(-m, 1, 1/2, l/(2x - 1)) 



m=0 



tt 2 (2x - 1) V2x-2 



>3C 



2Am \ p(-l/2,l/2-m) L - 



in- 



f vr (2x - l)y/2x - 2 m V 2x-l 



« 2A m [1 - l/(2s - l)]"" 1 d™ r r^-rn-l/2 
Y m\n (2x - l) 2 dz m L 



2=l/(2x-l) 



(48 



where P^ ,u (x) are the Jacobi polynomials. 

With this expression for X(x) the second integral in the third equality of Eq.(37) reads; 

] roo poo °° °° 

-/ dx dyX(x)\n\x -y\X(y) = ^ ^ A m A n L mn , (49) 

^ Jl Jl m=0n=0 
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where the matrix elements L mn = L nm are 



1 1 

2t ~ 2s 



(50) 



The third integral, with the aid of Eq.(45), reads 



- dx dyp(x) \n\x- y\p(y) = 2 r m r n R mn , 

* JO JO ™_n„-n 



(51) 



m=0 ri=0 



where 



-ftmn ■= - I dx I ay m x — y = . 

2 Jo io vx y/y 



(52) 



The first integral is, 



-1 rOC ^ ^ 

dx dyp(x) In |x - y\X(y) = V V r m A„C/ mn , 
JO Jl _n rr^ 



(53) 



m=0 ra=0 



where 



f^mn := / ds [ dt 
Jo Jo 



1 - s 



pM/2,l/2-n) (1 _ 2s) ^ 



1 /l 



2 Vs 



- (2t-l) 



(54) 



As an example we determine the leading term of — lnEp[J] of the \— ensemble. The 
density in this case is given by Eq.(23) with (5 — 2. We find the number of excluded 
charges in the interval (0, a), 



N(a)= dxa(x) ~ — y/a\nN, N -> oo. 



7T 



(55) 



In the limit N — > oo and a — > such that Y/alniV is finite, we find according to the 
screening theory that 



AF ~ s, iV — > oo, s := a [IniV] 2 finite, 



(56) 



which agrees with the results obtained from the un- folding transformation argument. 
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IV Solution via the continuum approach for a. — 1/2 

For a — 1/2 the solution of the integral equation reads 



1 b — x [ b dy y — a 1 



x r 

a J a 



n 2 (3 V x — a J a y — x\j b — y y/y 



b-x / . , x — a / b — a \\ . 7 . . 

k(p) + t — n ,p )), xe(a,b), (57) 



7i 2 pVb V x — a \ b — x \b — x 

where 

v( \ f 1 dy 2 b ~ a 

K(p) := / . p := — — 

Jo y/(l-y*)(l-p*y*) b 

is the complete elliptic function of the first kind, and 

f 1 dy 1 b- a 

U[n,p) := / r , n:~- 



lo 1-ny 2 yj(l-y*)(l-p*y2y ' 
is the complete elliptic function of the third kind. Note that the a(x) vanishes at x = b 
and is positive for x G (a, b). The upper band edge b is related to a and iV via the 
normalization condition; f% dxa(x) = N, as 

N=-^-y/E(K{p)-D{p)), (58) 

Tip 

where 



D(x) := - C dtJ- 4 o\ = -F(l/2,3/2,2,x 2 ) 

K ' 2 Jo V (l-t)(l-tx 2 ) 4 v 1 ' 1 ' ' ; 

is an elliptic function. The derivation of Eq. (58) is placed in Appendix A. In order to 
determine the free energy in the interval (a,b), we require the chemical potential A and 
the interaction energy | J^dxa(x). The chemical potential reads 

A = 2\fb - P [ b dxa(x) In \b - x\ = (3N (2 In 2 - In 6 - hip 2 ) +-Vbp 2 E(p), (59) 

where a derivation of the above is placed in appendix B. The interaction energy reads 

\ [ b dxu(x)a(x) = -\-{b - a) + -^-(6 - a) dt\D (VT^t)} 2 . (60) 

2 J a TT z p TT z p JO 1 v 7 J 
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The derivation of Eq.(60) is placed in appendix C. 



The free energy for N charges residing in (a, b) is recapitulated as follows: 




2NVbp 2 



E (p) + ^i b - a ] 1 + dtD 2 (VT^t) , 

(61) 



7T 



where 



N = 



iVb r 



E{p)--K{p) . 



(3ir I 



We now compute the change in the free energy, by expressing b as function of N and a in 
the limits iV — > oo and a — > 0, after some straightforward but tedious calculations (not 
reproduced here), we find 



which substantiates the result obtained from the screening theory. 
V Conclusion 

We have shown that applying the screening theory to the Coulomb fluid model of random 
matrix ensemble, gap formation distribution, which is the probability that there are no 
eigenvalue in a certain interval, J, of the spectrum is, in the limit of large interval 



Af[J] = J dxa(x) = number of eigenvalues excluded in J. 

This result shows that the determination of the leading term in —\nEp[J], requires only 
the integral of the local density cr(x), which enables us to make prediction on the gap 
formation probability of the a— ensemble at the hard edge. 



AF ~ — v/alniV 



Efs[J] ~ e 



-AT 2 




where 
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It is shown that for a > 1/2, the leading term in the asymptotic expansion is independent 
of a with the scaling variable, s oc Na. For a = 1/2, due to the logarithmic correction 
of the hard edge density the scaling variable, s oc [lniV] 2 a. However, for a < 1/2, the 
hard edge density acquires a dependence, and gives s oc N°a. Universality is violated for 



a < 1/2. Similar behaviour is found in the q— Hermite and q— Laguerre ensembles [pO 
The advantage of this approach is that it by-passes the requirement of the knowledge of 
the asymptotic expansion of non-classical orthogonal polynomials, especially those asso- 
ciated with the indeterminate classical moment problems. 

As a side product of the screening theory, we make prediction on the level spacing distri- 
bution in the soft edge of Hermitian random matrix models that arises in double-scaling 
limit of the two-dimensional quantum gravity. For example, taking u(x) = N(x + j^£ 4 ), 



it can be shown that [KB] tuning g to g c (< 0), the critical point density reads, 

V ; 24vr ; ' 

where now J^gdxa(x) = 1. Note that in this formulation, 

F \ rb rb rb ( 9 l 

— = -2j_ b dx J_ dycr(x) \n\x-y\a(y) + J^dxa{x) f x + —x l 
Using our theory, a simple calculation gives, 

-\nE 2 (s) ~ N' 



V8 



dxaJx) 



~s 5 



where the scaling variable is s := A^ 2 / 5 (v^8 — a), in the limit N — > oo, a — >• y/8 such that s 
is finite. The detail exposition for this and other multi-critical point models can be found 



m 



29| . Work is underway to use this approach to determine, Ep(s), of the q— Laguerre 
ensembles. 
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Appendix A The normalization integral 

By definition, 

f b , /x I b - a r 1 1 I 1-y r 1 dx I x 

N = L dxa(x) = Jo dy yw=rt)h ^VT^ 

1 b — a f 1 , / 1 — v 1 b — a , , , . . . , 9n 

= -^^T / d V\~h —\=^-HT B V2,3/2 F 1/2, 1/2,2, 

= A vV[^(p)-^(p)], (ai) 

where we have made of 

F(l/2, 1/2, 2,p 2 ) = - [K(p) - Dip)] . (A2) 
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Appendix B Evaluation of the chemical potential 



We have 



where 



A = 2Vb - (3N\n(b - a) + I(a,b), 

I(a,b) := —(3(b — a) / dx\nxa(b — x(b — a)) 
Jo 

2\fb 2 r 1 dy I 1 - y r 1 dx I x~ 
7r 2 ^ Jo 2^/y\ 1 — p 2 y Jo x — y y 1 — x 

2^2(1 I l-y 2 d f 1 dx x"- 1 



(B.l) 



lnx 



7r 2 7o V i — v 2 v 2 J° 



X 



y 2 yJl — X 



j/=3/2 



We now focus on the derivative with respect to u, 



u=3/2 



(B.2) 



|-5(-l/2, u)F(3/2 -u,l, 3/2, 1 - y 2 ) 



d 



u=3/2 



F(0, 1, 3/2, 1 - y 2 ) +B(-l/2, 3/2)-^F(3/2 - u, 1, 3/2, 1 - y 2 ) 

u=3/2 OV 



= -2tt(1 - In 2) + ^F(z, 1, 3/2, 1 - y 2 ) 



z=0 



u=3/2 

B.3) 



Note that we can make use of 



d T{z + a) T{z + a) 



{z + a) -ip(z + b)], 



dzT(z + b) T(z + b) 
to compute the derivative of the Beta function with respect to its argument. We now 
only have to compute the derivative of the Hypergeometric function with respect to its 
first parameter. As z + 1 — 3/2 < 0, for z sufficiently close to 0, we may employ the 
power series representation of F(z, 1, 3/2, 1 — y 2 ) which is valid in < |1 — y 2 \ < 1 and 
differentiate term by term resulting in 



^F(z, 1,3/2, l-y 2 ) 



2\n 



= E%3/2)(1-^) 

2 = n = 1 
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/•i dt ^ r 
■/o tVl - i„=i L 

( _ 2) f 1 dt 

= [ V } Jo VT^t[l-t(l-yi)] 

= (l-y 2 )5(l,3/2)F(l,l, 5/2, 1 - y 2 ) 



= 2 



2?/ 



From Eqs.(B.4),(B.3) and (B.2) we find 



sin^/T^). 



(S.4) 



I(a,b) 



2Vb 



7T^ 



(5.5) 



The integrals in Eq.(B.5) are recognized as the standard elliptic integrals (K and D), and 
we find, 



I(a, b) = 2Vb 



7T 



7T 



7T 



(5.6) 



Observe that D(p) can be expressed in terms of N and K(jp) [see Eq. (A.l)]. Further 



more using 



p 2 D(p) = K(p)-E(p), 



where 



2„2 



we have, 



and therefore 



/(a, 6) = 2/3iV In 2 - 2Vb + —p 2 E(p) 



A\/b 

A = /37V(21n2 - In 6 - hip 2 ) + —p 2 E{p). 
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Appendix C Evaluation of the interaction energy 

In this appendix we determine the interaction energy for a > 0, thus, 



rb 

/ u(x)a(x) — h(a, b, a) + h{a, b, a), 

J a 



l r b 
2 



where 



h(a, b, a) := 



b^ib-a 
27r 2 a(3 



r d x rj^ MLzi f, 

Jo Jo x — y y — x 



p 2 y 



2a- 1 



and 



r , , \ & 2Q_1 (&-a) f 1 , f 1 dy 
i2{a,b,a) := — ^-^ — - — / ax 



2n 2 af3 Jo Jo x — y 



We first compute h^a, b, a), 



J 2 (a, b, a) 



a 2a - 1 (6-a) r 1 
2na(3 Jo 



dy(l-p y) 



y 2(5 2a 



F(l 



and find 



7 2 (a,6,l/2) = 



b — a 



We now compute ii(a, 6, a). Observe that, 



\i-P 2 y) 


1 


y-x 




J p 2 (x-y)" 
1-p 2 ?/ _ 


- 1 


2/ - x 



-L 



i- p 2y de d 
x — ydt 



[l-e(x-y)] c 



a-l 



With the change of variable 



1 + ey 
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t, 



we get 



\i-p 2 yj 



y - x 

Using Eq. (C.5) in Eq. (C.2), 



rv 2 

a 

Jo 



dt 

(1 -ty) 1+a (l -txy- a ' 



(C.5) 



h(a, b, a) = 



b 2a ~\b 



x f 1 l — yfp 



dy 



dt 



2tt 2 [3 Jo ]] 1- x Jo ^ y Jo (1 - ty) 1+a (l - tx) 1 -* 

b 2a -\b-a) o/oA /-f 2 JiJ71/1 o/o f 1 ^{i - p 2 y) 2a ~ l l^y 



For a = 1/2, Eq. (C.6) becomes 

b — a 



5(1/2,3/2) T dtF(l-a,3/2,2,t) I dy^— 
Jo Jo 1 



ty 



l+a 



(C6) 



h(a, b, 1/2) = -_^£ 2 (l/2, 3/2) jT d*F(l/2, 3/2, 2, t)F(3/2, 1/2, 2, f) 



b- a fP 2 



8/3 



(C.7) 



where we have used F(3/2, 1/2, 2, t) = F(l/2, 3/2, 2, t) = (v^) , to get Eq. (C.7). 
From Eqs. (C.4) and (C.7), the interaction energy reads 



/ 1 (o,6,l/2)+/ 2 (o,6,l/2) = ^U 



A^dt[D(Vi)] 2 ). (C.8) 



Note that 



£ dt[D (Vi)f = £ dt[D (Vi)f - dt[D (Vi)] 2 , (C.9) 



2 r a / b 



since p 2 = 1 — a/b. With the aid of 



[dt[D(Vt) 



2 7T 
= T ' 



and Eq. (C.9), we arrive at Eq. (60) in the text. 
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